Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer.

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Whole-genome analysis of human papillomavirus 67 isolated from Japanese women with cervical lesions.

BACKGROUND: Human papillomavirus (HPV) type 67 is phylogenetically classified into Alphapapillomavirus species 9 (alpha-9) together with other carcinogenic types (HPV16, 31, 33, 35, 52 and 58), but is the only alpha-9 type defined as possibly carcinogenic. This study aimed to determine whole-genome sequences of HPV67 isolated from Japanese women with cervical cancer or cervical intraepithelial neoplasia (CIN) to better understand the genetic basis of the oncogenic potential of HPV67. METHODS: Total cellular DNA isolated from cervical exfoliated cells that were single positive for HPV67 (invasive cervical cancer, n = 2; CIN3, n = 6; CIN 2, n = 1; CIN1, n = 2; the normal cervix, n = 1) was subjected to PCR to amplify HPV67 DNA, followed by next generation sequencing to determine the complete viral genome sequences. Variant lineages/sublineages were assigned through viral whole-genome phylogenetic analysis. The transcriptional activity of the virus early promoter was assessed by luciferase reporter assays in cervical cancer cell lines. RESULTS: Phylogenetic analyses of HPV67 genomes from Japan (n = 12) revealed that 11 belonged to lineage A (sublineage A1, n = 9; sublineage A2, n = 2) and one belonged to lineage B. All cancer cases contained sublineage A1 variants, and one of these contained an in-frame deletion in the E2 gene. The long control region of the HPV67 genome did not induce transcription from the virus early promoter in HeLa cells, but showed weak transcriptional activity in CaSki cells. CONCLUSIONS: The single detection of HPV67 in cervical cancer and precancer specimens strongly suggests the carcinogenic potential of this rare genotype. The phylogenetic analysis indicates a predominance of lineage A variants among HPV67 infections in Japan. Since only 23 complete genome sequences of HPV67 had been obtained until now, the newly determined genome sequences in this study are expected to contribute to further functional and evolutionary studies on the genetic diversity of HPV67.

Ancient Evolutionary History of Human Papillomavirus Type 16, 18 and 58 Variants Prevalent Exclusively in Japan.

Human papillomavirus (HPV) is a sexually transmitted virus with an approximately 8-kilo base DNA genome, which establishes long-term persistent infection in anogenital tissues. High levels of genetic variations, including viral genotypes and intra-type variants, have been described for HPV genomes, together with geographical differences in the distribution of genotypes and variants. Here, by employing a maximum likelihood method, we performed phylogenetic analyses of the complete genome sequences of HPV16, HPV18 and HPV58 available from GenBank (n = 627, 146 and 157, respectively). We found several characteristic clusters that exclusively contain HPV genomes from Japan: two for HPV16 (sublineages A4 and A5), one for HPV18 (sublineage A1) and two for HPV58 (sublineages A1 and A2). Bayesian phylogenetic analyses of concatenated viral gene sequences showed that divergence of the most recent common ancestor of these Japan-specific clades was estimated to have occurred ~98,000 years before present (YBP) for HPV16 A4, ~39,000 YBP for HPV16 A5, ~38,000 YBP for HPV18 A1, ~26,000 for HPV58 A1 and ~25,000 YBP for HPV58 A2. This estimated timeframe for the divergence of the Japan-specific clades suggests that the introduction of these HPV variants into the Japanese archipelago dates back to at least ~25,000 YBP and provides a scenario of virus co-migration with ancestral Japanese populations from continental Asia during the Upper Paleolithic period.